Fluid Dynamics of Relativistic Quantum Dust 



Hans-Thomas Elze 

I Universidade Federal do Rio de Janeiro, Instituto de Fisica 

§ ■ Caixa Postal 68.528, 21945-970 Rio de Janeiro, RJ, Brazil 

(N ■ 

j|H ' December 2001 

■ The microscopic transport equations for free fields are solved using the Schwinger function. Thus, 
I for general initial conditions, the evolution of the energy-momentum tensor is obtained, incor- 

■ porating the quantum effects exactly. The result for relativistic fermions differs from classical 
I hydrodynamics, which is illustrated for Landau and Bjorken type initial conditions in this model 

of exploding primordial matter. Free fermions behave like classical dust concerning hydrodynamic 
observables. However, quantum effects which are present in the initial state are preserved. 



Abstract 



o 

■ PACS numbers: 



Often the complicated time dependent dynamics of quantum many-body systems or fields is approx- 
imated by a perfect fluid model. Since the seminal work by Fermi and Landau this approach has 



I been applied successfully, in order to study global features, such as multiplicity distributions and 

• apparently thermal transverse momentum spectra of produced particles, in high-energy collisions of 

strongly interacting matter [||, |2|, ^, |^. Similarly, the hydrodynamic approximation is often invoked 
^ \ in astrophysical applications and cosmological studies of the early universe |^ 

\^ • Recently it has been shown that a free scalar field indeed behaves like a perfect fluid in the semiclas- 

^ sical (WKB) regime More generally, the mechanisms of quantum decoherence and thermalization 

in such systems which can be described hydrodynamically, i.e. the emergence of classical deterministic 
evolution from an underlying quantum field theory, are of fundamental interest |^, ^, ^, 10|. 



The limitations of the fluid picture, however, have rarely been explored in the microscopic or high 
energy density domain. Difficulties reside in the derivation of consistent transport equations and in the 
amount of computation required to find realistic solutions; see Refs. [pT|], for example, for a review and 
recent progress concerning selfinteracting scalar particles and the quark-gluon plasma, respectively. 
More understanding of related hydrodynamic behavior, if any, seems highly desirable. 

Presently, we study the relation between relativistic hydrodynamics and the full quantum evolution 
of a free matter field. In the absence of interactions, decoherence or thermalization may be present in 
the initial state, corresponding to an impure density matrix, but is followed by unitary evolution. We 
consider this as a "quantum dust" model of the expansion of matter originating from a high energy 
density preparation phase, which the Landau and Bjorken models describe classically ||2|, ^. 

Our approach is independent of the nature of the field, as long as it obeys a standard wave 
equation. To be definite, we choose to work with Dirac fermions and comment about neutrinos 
later. We introduce the spinor Wigner function, i.e., a (4x4)-matrix depending on space-time and 
four-momentum coordinates: 

W^p{x;p) ^ J ^e-'P-y{:M^ + y/2)Mx-y/2) :) , (1) 
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where the expectation value refers to the (mixed) state of the system; without interactions, the vacuum 
plays only a passive role and, therefore, is eliminated by normal-ordering the field operators. 

All observables can be expressed in terms of the Wigner function here. In particular, the (unsym- 
metrized) energy-momentum tensor: 



{■.T^,{x):)=i{:^{xh^d,^{x):)=tT^^J d''pp,W{x;p) , (2) 

where d= — d) and with a trace over spinor indices (conventions as in Furthermore, the 

dynamics of W reduces to the usual phase space description in the classical limit |T2| . 

Propagation of the free fields entering in Eq. (||) from one time-like hypersurface to another is 
described by the Schwinger function. It is the solution of the homogeneous Dirac equation, [ij ■ 
dx — m\S{x,x') = 0, for the initial condition S{x,x',x^ = x'^) = —ij^6^{x — x'). Thus, ip{x) = 
i J d'^x' S{x, x')^'^ip{x'), and similarly for the adjoint. An explicit form is: 

iS{x,x') = iS{x -x'^A) = {ij -d^ + m) j -^0^^ - e-^^^-^) , (3) 



where k± = (iwfc, k) and uJk = {k^ + m^)^/^. 

Making use of Eqs. (^ and (|3|), we relate the Wigner function at different times, t 



,0 uo. 



W{x-p) = J -^e-^^--5^{pt)5^{pl) j d?x'e^^--' j dp'^K{pt)fW(x';p')^^K{pl) , (4) 

where = p it |, (5^(g) = ±(5((?^ — vn?) (for /\q^\ = ±1), A(g) = ^ ■ q + m, and p'^ = {p'^,p)- 

The Eq. (^) implies that the Wigner function obeys a generalized mass-shell constraint and a proper 
free-streaming transport equation: 

\p^-m^-—dl]W{x;p) = , (5) 

p-dx Wix;p) = , (6) 

separately for each matrix element. The reinserted fi indicates the important quantum term in the 
equations, which otherwise have the familiar classical appearance. 

Thus Eq. (^) presents an integral solution of the microscopic transport equations for a given initial 
Wigner function. Furthermore, a semiclassical approximation of the Schwinger function may be used 
to generate an integral solution of the corresponding classical transport problem. 

Next, we decompose the Wigner function with respect to the standard basis of the Clifford algebra, 

■^(T^'^S^yj i.e., in terms of scalar, pseudoscalar, vector, axial vector, 
and antisymmetric tensor components. The functions, T = |tr W = —\i tr '^^W , = \ti "j^iW , 
A/j, = |tr 757^VF , and 5^,^ = ^tr a^^W , which represent physical current densities, are real, due to 

= 7°VF7° Q. They individually obey Eqs. (|) and (|). 

We assume "P = = Afj,, i.e., we consider a spin saturated system for simplicity Then, 
using the 'transport equation' which follows directly from the Dirac equation applied to W ^ [7 ■ (p + 
^dx) — m]W{x]p) = 0, and decomposing it accordingly, the following additional relations among the 
remaining densities are obtained: 

v^{x-^p) = "^n^-.p) , (7) 

s^'^ix-p) = ^(j^'d^-p'd^x)Hx;p) ■ (8) 
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Note that S^'^ is intrinsically by one order in h smaller than the other two densities. 

We conclude that presently the dynamics of the system is represented completely by the scalar 
phase space density T. Using Eqs. (^) and (j^, we obtain in particular: 

(: T'^'ix) :)=4mJ d^p ?^T{x;p) , (9) 

which is symmetric and conserved, d^T^^{x) = 0, on account of Eq. (^). Furthermore, this implies 
the 'equation of state': 

(: rOO(x) :) - ^(: :) = Am j d^p , (10) 

i=i •' 

which relates energy density and pressure(s). However, applying Eq. (^), we find that this relationship 
evolves in a wavelike manner, driven by off-shell contributions to the evolving T: 



d^{: Tjiix) :) = 16m j d^ {p - m')T{x;p) . (11) 
This differs from classical hydrodynamics with a fixed functional form of the equation of state. Eqs. (0)- 



(|ll| ) hold independently of the initial state, of course, if it evolves without further interaction. 

Making use of Eq. (^ in Eq. (^) , we now calculate the energy-momentum tensor at any time in terms 
of the initial scalar density. Employing the decomposition of the Wigner function and commutation 
and trace relations for the 7 matrices, as well as Eqs. (^)-(|8|), we obtain: 



(: r-(x) :) = 8m / . . . / ^ / + ^ + P') - ^ {ip'ky - p^P 



+p'p"k^ + ^ {k')Vfj Hx';p,p") , (12) 

where J ■ ■ ■ J = (27r)~^ / d^ / d^/c e~*'^'^J^(p^)5^(p^) / d^x'e*^'^' / dp'^; we also made use of partial 
integrations and the (5-function constraints. The three terms on the right-hand side stem from the 
scalar, vector, and antisymmetric tensor components of the initial Wigner function, respectively. 

If the initial distribution is an isotropic function of the three-momentum, then T^"^ is diagonal at 
all times, implying that the absence of flow in the initial state will be preserved. 

Indeed, we expect the (non-)flow features of the initial distribution to be preserved during the 
evolution, due to the absence of interactions. Kinetic energy from microscopic particle degrees of 
freedom will not be converted into collective motion. An interesting question is, how the classical hy- 
drodynamic acceleration of fluid cells due to pressure gradients arises in our present model after coarse 
graining I, I, I, |10|. We do not pursue this at present. Recalling earlier work on the hydrodynamic 



representation of quantum mechanics, e.g. Refs. [0], and recently deduced classical fluid behavior of 
quantum fields in WKB approximation [Q, we study the full quantum effects here. 

We consider the exact evolution of T^'^, assuming a particle-antiparticle symmetric initial state. 
This is believed to hold, for example, close to midrapidity in the center-of-mass frame of central 
high-energy collisions It implies that the initial J- is an even function of the energy variable, 

J-{x';p,p'^) = J-{x';p, —p'^). While Eq. ([T^) allows general initial conditions, we follow the implicit 
on-shell assumption in classical hydrodynamic models: 

J^{x';p,p'^) = (27r)-3m S{p'^ - m^) (e{p'^)F{x';p,p'°) + e{-p'^)F{x';p, -p'°)) . (13) 

Fermion blackbody radiation is described by F{x';p,p'^) = f{p'^/T{x')), where T denotes the local 
temperature, and with f{s) = (c^ -|- 1)^^; this is easily illustrated with the help of Eqs. (P) and (p^). 
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Implementing Eq. (^), we obtain the simpler result: 

+ uj-)'^ — k^) cos[(w+ — uJ-)t\ — ~ '^-)^ ~ cos[(w+ + (J-'-)*]} , 

where t = — x"^, ujp = {p^ + m^)^/^, l<j± = ((p it A;/2)^ + m^)-*^/^; furthermore, p'^ = ± 
with "+" when multiplying the first and "— " when multiplying the second term of the difference, 
respectively. Depending on geometry and initial state, further integrations can be done analytically. 

Consider a (l+l)-dimensional system for illustration, assuming that the particles are approximately 
massless, i.e. ojp ~ and that F is even inp (no flow). Specializing to a Landau type initial condition, 
the distribution is prepared on a fixed timelike hypersurface at t = [Q. We find the ultrarelativistic 
equation of state for the only nonvanishing components of T^'^, e = T^^ = T^^ = P {d = 1 + 1), which 
are calculated as a momentum integral following Eq. (p!4[): 

T^\x,t) = 2l^\p\[F{x-t;\p\)+Fix + t;\p\)) (15) 

t,t = 0) +T™(x + t,t = 0)) , (16) 



1 (T"(. 



i.e., a superposition of wavelike propagating momentum contributions in accordance with Eq. (|l l]). 

Similarly, a Bjorken type initial condition can be specified on a surface of constant proper time 
Q. A transformation of Eq. (|l^) to space-time rapidity and proper time coordinates yields: 

r°0(y,T) = 2 J ^|p|(F(-Toe-J^/2+ln.Ao. + i7(^„e?//2+lnr/ro. IpD) ^ (17) 

since x = rsinhy/2 and t = TCOshy/2 (r > tq > 0). 

Our results for free fermions show the free-streaming behavior of classical dust, associated with 
the independent propagation and linear superposition of the momentum contributions to the scalar 
component F of the Wigner function here. In particular, the shape function of each mode is preserved 
and translated lightlike (with dispersion for massive particles). Due to the assumed momentum sym- 
metry, the initial distribution will separate into two components after a finite time, travelling into the 
forward and backward direction, respectively, with a corresponding dilution at the center. 

We recall that T^"^ being diagonal implies the absence of ideal hydrodynamic flow, given e = 



{d— 1)P- This does not depend on whether the initial state is on- or off-shell, see Eq.(12). Therefore, 
any hydrodynamic behavior must be the effect of a pecularity of the semiclassical limit [^], of coarse 
graining [^, |l^, or of interactions [^], or a combination of these. 

Despite the apparently classical evolution, however, all initial state quantum effects are incor- 
porated and preserved. If the initial dimensionless distribution F has a dependence on products of 
momentum and space-time variables, which is characteristic for matter waves, such terms invoke a fac- 
tor 1/h. Similarly, if it is thermal (T) but includes the finite size (L) shell effects or global constraints, 
then there are quantum corrections involving LT/h (ks = c = 1) |l^. They have not been included in 
semiclassical tranport or classical hydrodynamic models of high-energy (nuclear) collisions, but may 
be large. Here the quantum dust model provides a testing ground to assess the importance of these 
quantum effects. 

Our approach based on the Schwinger function reduces the solution of the free quantum transport 
problem to quadratures - analytical results in three dimensions can be obtained and will be discussed 
elsewhere. It may lead to an efficient way of treating interacting particles, when a perturbative 
expansion is meaningful. 
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